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Abstract 

We introduce and discuss spatially homogeneous Maxwell-type models of the 
nonlinear Boltzmann equation undergoing binary collisions with a random com- 
ponent. The random contribution to collisions is such that the usual coUisional 
invariants of mass, momentum and energy do not hold pointwise, even if they all 
hold in the mean. Under this assumption it is shown that, while the Boltzmann 
equation has the usual conserved quantities, it possesses a steady state with power- 
like tails for certain random variables. A similar situation occurs in kinetic models 
of economy recently considered by two of the authors |24] , which arc conservative in 
the mean but possess a steady distribution with Parcto tails. The convolution-like 
gain operator is subsequently shown to have good contraction/expansion proper- 
ties with respect to different metrics in the set of probability measures. Existence 
and regularity of isotropic stationary states is shown directly by constructing con- 
verging iteration sequences as done in [5]. Uniqueness, asymptotic stability and 
estimates of overpopulated high energy tails of the steady profile are derived from 
the basic property of contraction/expansion of metrics. For general initial condi- 
tions the solutions of the Boltzmann equation are then proved to converge with 
computable rate as t — > cxd to the steady solution in these distances, which metri- 
cizes the weak convergence of measures. These results show that power-like tails 
in Maxwell models are obtained when the point-wise conservation of momentum 
and/or energy holds only globally. 

1 Introduction 

In this paper, we introduce and discuss the possibihty to obtain steady solutions with 
power-like tails starting from conservative molecular systems described by the Boltz- 
mann equation with Maxwell-type collision kernels. The starting point of our model is 
to consider binary collisions that result in a linear combination of an inelastic collision 
and a random contribution. As we shall see, the random addition to the post-collision 
velocities can only increase the mean of the collisional energy, and, among other things. 
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it gives the possibility to construct a binary collision that preserves (in the mean) mass, 
momentum and energy. Our model is closely related to a kinetic model for economics 
introduced by Pareschi and two of the present authors [18]. There, the random contri- 
bution to the collision (trade) was introduced to take into account the returns of the 
market. 

Inelastic Maxwell models were introduced by Bobylev, Gamba and one of the au- 
thors in 2000 [6]; see also [23] for the one dimensional case. Maybe the most interesting 
result (absent in the elastic case) is the existence of self-similar solutions in the ho- 
mogeneous cooling problem and the non-Maxwellian behavior of these solutions, which 
displays power-like decay for large velocities. It was conjectured in [20] and later proved 
in [U \TT\ [2] that such solutions represent intermediate asymptotics for a wide class of 
initial data. Other results concerned with self-similar solutions in the theory of the 
classical (elastic) Boltzmann equation for Maxwell molecules were also recently pub- 
lished in [318]. In light of these results, it looks clear that in many aspects both elastic 
and inelastic Maxwell models must be studied from a unified point of view. As ob- 
served by Bobylev and Gamba in [9l [10] , an interesting question arises in connection 
with power-like tails for high velocities. Is it possible to observe a similar effect, i.e., 
an appearance of power-like tails from initial data with exponential tails, in a system 
of particles interacting according to laws of classical mechanics without energy loss? In 
|12j Bobylev and Gamba gave a partial answer to this question by showing that, under 
a certain limiting procedure, such behavior can in principle be observed if one considers 
a mixture of classical Maxwell gases. More precisely, self-similar solutions converging 
towards maxwellian equilibrium were proved to have power-like tails once normalized 
by the equilibrium. 

In this paper, we will try to elucidate the same question, starting from a somewhat 
different point of view. Our starting point will be a suitable modification to the ho- 
mogeneous Boltzmann equation for the inelastic Maxwell molecules introduced in Ref. 
[6], in such a way that the usual conservations of mass, momentum and energy in the 
binary collisions still continue to hold in the mean sense. The scaled-in-time inelastic 
Boltzmann equation introduced in [6] reads 

% = Qe{f,f). (1) 

Here, f{v, t) is the density for the velocity space distribution of the molecules at time t, 
while Qe{f,f) is the inelastic Boltzmann collision operator, which contains the effects 
of binary collisions of grains. As usual in this context, the collision operator Qeif,f) 
is more easily treated if expressed in weak form. This corresponds to writing, for every 
suitable test function ip, 

{v,Qe{f,f)) = ^ f [ [ f{v)f{w)L{v*)-ipiv)]dvdwda. (2) 
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In ([2]), V* is the outgoing velocity corresponding to a particle in the collision defined 
by the incoming velocities v,w and the angular parameter a ^ S'^: 

* 1/ N , 1 + e, , 

V =-(v + w]-\ (v — w] -\ \v — w\a, 

(3) 

* 1/ N , 1 + e, , 

^' = 2 '•^ ~'~ 4 — '-^ — wj — \v — w\a . 
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The parameter < e < 1 represents the restitution coefficient. 

In the model we consider, this restitution coefficient will be chosen as a random 
variable that can be interpreted from a physical point of view as the stochasticity in 
the microscopic process of collision due to the randomness of the grains' geometry and 
the mechanical properties of the medium. We will show in the next section that this 
random behavior in restoring energy leads to a precise form of the energy gain term 
that differs from the usually chosen diffusion term, the so-called "thermal bath". This 
new form of "thermal bath" is thus related to the process generated by the randomness 
of the granular media. We prove that this particular thermal bath yields equilibrium 
states with power law tails. 

Such over-populated tails in distributions at equilibrium arise in other contexts. 
We shall present similar results on the large time behavior of collisional kinetic the- 
ory applied to economic modelling. In this framework, the kinetic variable represents 
the wealth of agents and the collision operator describes the evolution of the wealth 
distribution through exchanges. We refer to [18\ [2l] and references therein for a math- 
ematical presentation of these models closely related to so called " econo-physics" . In 
such models, the equations between pre- and post-collisional values involve some ran- 
domness that is related to the stochasticity of the market that provides random returns. 

In the remainder of this paper, we will study in detail the large time behavior of the 
solution of the Boltzmann equation involving such a stochastic process. We show that 
the validity (at a macroscopic level) of the classical collision invariants is enough to 
guarantee convergence towards a steady profile, but not enough to reach a Maxwellian- 
like profile. In fact, we will show that there is a class of random perturbations of the 
coefficient of restitution such that the steady state possesses power-like tails. 

A crucial role in our analysis is played by the weak norm convergence, which is 
obtained by further pushing the development of a method first used in [21] to control 
the exponential convergence of Maxwellian molecules in certain weak norms. This 
will be done by using the fact that the nonlinear operator in the Boltzmann equation 
(see ([9])) can be expressed in Fourier variables in a simple closed form using Bobylev's 
identity [5] . Estimates of the evolution of the Wasserstein distance \32\ [29l [30] between 
solutions will be presented for the economic and the inelastic model since they give 
complementary information with respect to the results in [23]. Concerning this second 
aspect, we will take advantage of the recent analysis of Bolley and Carrillo |14[ll7j of the 
inelastic Boltzmann equation for Maxwell molecules. From this analysis, we will obtain 
the uniqueness and asymptotic stability of stationary states for this model. Finally, 
the appearance of power-like tails for the asymptotically stable stationary states will 
be discussed for both models, giving explicit examples of random variables producing 
this behavior. 

The paper is organized as follows: in section [2] we detail the collisional models for 
both granular media and economy applications including random coefficients in the 
relations between pre- and post-collisional variables. In section [3l we recall the main 
properties of probability metrics. In section HJ we investigate large time behavior of the 
solution of the kinetic economy model and section [5] is devoted to large time behavior 
of stochastic granular media. 

Lastly, let us summarize the two main results of this paper : first, we give some 
insight into conditions for a collision operator to lead to power-law tails (conservatism 
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in mean being some kind of necessary condition); second, we propose a new form for 
the thermal bath with a physicahy relevant origin (the restitution coefficient taking 
into account the randomness of granular media). 

2 Modelling issues and diffusion approximation 

Let us present the proposed stochastic granular model (with a random restitution 
coefficient) and its diffusion limit and then recall briefly the similar analysis for the 
economy model following |18| . 

2.1 Stochastic granular media 

Considering the weak formulation ([2]), easy computations show that {ip{v), Qdf, /)) = 
whenever ip{v) = 1 and ip{v) = v, while {^iv),Qeif,f)) < if ip{v) = . This 
corresponds to conservation of mass and momentum, and, respectively, to loss of energy 
for the solution to equation ([1]). For this reason, if we fix the initial data to be a centered 
probability density function, the solution will remain centered at any subsequent time 
t > 0. The loss of energy in a single collision with a constant restitution coefficient e is 
given by 

1 _ e2 

+ \w'^ = |f p + |tt;p ^ — (|f — u)p — \v — w\{v — w) ■ . (4) 

The previous formula is the key to our modification of the collisions. Let us replace the 
constant coefficient of restitution e with a stochastic coefficient of restitution e, such 
that for a given random variable r/ 

e = e + ?7, with (r/) = and {rf') = (5^ . (5) 

In dS]) and in the rest of the paper, (•) denotes the mathematical expectation of the real- 
valued random variable ry, i.e., integration against a measure [i. For obvious physical 
reasons, the random variable r/ has to be chosen to satisfy 77 > — e, in order to guarantee 
that the (random) coefficient of restitution e > 0. Using e instead of e in ([3]) gives that 
the momentum is conserved in average for a suitable choice of the variance. In fact, 
since 

+ = + {wl"^ — — (jv — — \ v — w\{v — w) ■ a) , (6) 

by choosing the variance we obtain 

(|^;f + |«;'|2) = |^;|2 + |^|2. (7) 

We will call a collision process (or equivalently a random cross section) satisfying ([7]) 
conservative in the mean. Let us remark that condition d?]) cannot be satisfied if e takes 
only values less than 1, since in that case remains also less than 1 and so does its 
average < >= + < 1. The main idea behind this is that particles can even 
gain energy in collisions even though the total energy is conserved in the mean. 

From the physical point of view, this assumption of energy-gain particle collisions 
may seem strange. We will show in the sequel that this energy input can be interpreted 
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as a sort of thermal bath. Particles are immersed in a medium that produces this 
random change in the strength of their relative velocity. We will argue, based on 
a derivation of a Fokker-Planck approximation, that this random component in the 
collision operator can be approximated by a second-order differential operator whose 
diffusion matrix depends on the second moments of the solution / itself and the random 
variable rj (see |18l [25] for a similar approach in one dimension) . 

This idea allows us to consider a new class of Maxwell-type models, from now on 
called conservative in the mean, which are obtained from (post-collision) velocities given 
by 

, 1, , 1 — e, , 1 + e. . 
V =-(v + w]-\ ■ — (v — w) -\ : — \v — w\a, 

(8) 

,1 1-e 
w =-(t' + 'wj ^ — [y — w) ^ — \v — w\a . 

where e is the random coefficient of restitution defined in ([5]), and = 1 — e^. The 
corresponding Boltzmann equation reads 

^=QeifJ) = {QeifJ)) , (9) 

and its corresponding weak form is 

iv>,Qe{f,f)) = ^( [ I [ fiv)fiw)\ipiv')-^iv)\dvdwda). (10) 
4vr \ Jk3 Jir3 Js2 I J / 

In view of our choice of the random contribution to the coefficient of restitution, we 
now have {^{v),Qe{f,f)) = whenever ip{v) = l,v, that is, the classical collision 
invariants of the elastic Boltzmann equation. 



2.2 Formal diffusive asymptotics 

Before entering into the study of the large-time behavior of the Boltzmann equation 
([9]), we shall present here some formal arguments that hopefully clarify the action of 
the random restitution coefficient in the collision mechanism ([8|). 

To this end, following the same method as in [2^, letting (v',w') denote the post- 
collision velocities ([8]) in our random collision with {v*,w*) as post collision velocities 
defined by the classic inelastic collision ([3]), we can split the velocities into their deter- 
ministic and random parts 

v' = V* + r]A{u,a) , w' = w* — r]A{u,o-), (11) 

where we let u = v — w and 

^{u,(^) = ^ {\u\o- - u) . 

Let us consider a Taylor expansion of (p{v') around f{v*) up to second order in r/. 
Thanks to (jlip we get 

^(y') = ^(^*) + ,j (V(^(^*) • A(n, a)) + E ^^^^^^ + " " " ^^^^ 
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Thus, taking the mean of the expansion (jl2p . and using the property {vj) = 0, we get 



(13) 



Truncating the expansion (|13p after the second-order term and inserting (jlSp into (jlOp , 
we conclude 



, Qe(/, /)) ^ (V^ , Qe(/, /)) + {V, Deif, /)) 



(14) 



(V', Qe(/,/)) + 



Svr 



JiRS dv*dv* 



While the first term in (jl4p Qe{f, f) is the classical inelastic Boltzmann collision oper- 
ator, the second term De{f, f) needs to be further analyzed. 

Denoting by {*v, *w) the pre-collision velocities in the inelastic collision, and taking 
into account the fact that the Jacobian of the transformation d*vd*w into dv dw for a 
constant restitution coefficient is equal to e~^, one obtains 



{if, D,{f,f)) 



Svr 



Stt 



52 ^ dv*dv* 



E 



AiAjf{v)f{w)dv dw da 



M3 Js2 e f-^ dvidvj 



ip{v) 



E 



9^V3(i 



dvidvj 



S2 e 



'Ai*Ajf{*v)f{*w)dwda 



dv 



Stt ^ dvidvj J ^3 Js^ e 



f{*v)f{*w)dwda 



dv. (15) 



This shows that, at least for small inelasticity, the random part of the collision corre- 
sponds to a correction given by the nonlinear diffusion operator De{f, f){v), where 

DeUJ){v) = ^Y.irjr I I -*\*^jfCv)fCw)dwda. (16) 
Svr ^ dvidvj J^s J 32 e 



Different expressions of the operator (I16p can be recovered owing to the definition of 
A. For the purposes of the present paper, however, we simply remark that, choosing 
the test function ip{v) = jfp, direct computations show that the correction D(,{f,f) is 
such that 



i\v\\D,if,f)) 



647r 



52 



|M|cr — n| f{v)f{w)dvdwda 



\v\'^f{v) dv 



V f{v) dv 



(17) 



This reveals the fundamental fact that the diffusion operator produces a growth of 
the second moment proportional to the second moment itself. This action is clearly 



6 



different from the action of a linear diffusion operator (a tliermal bath), which induces 
a growth of the second moment proportional to the mass. This supports the fact that 
the Boltzmann equation ([9]) can produce fat tails. 



2.3 Simple economy market modelling 

In one dimension of the "velocity" variable, a similar construction leads to kinetic 
models for wealth redistribution I24j . In this case, the variable v E R-|_ represents 
the wealth of the agents, binary collisions are trades between agents, and the (eventual) 
power-like tails of the steady distribution of wealth are known in the pertinent literature 
as Pareto tails. Due to the fact that the variable is in M_|_, the possible conserved 
quantities reduce to mass and momentum. In [18] the collision mechanism is given by 

v' = {1 - \)v + Xw + -qv; w' = Xv + {1 - \)w + 'n*w (18) 

where < A < 1 represents the constant saving rate and t] and r]* are random variables 
with law given by a measure fi{s) of zero mean, variance and support in [—A, +oo). 
In this way, for all realizations of the random variable we have rj > —X and wealths 
after trading are well defined i.e., remain nonnegative. This is the so-called no debt 
condition. In this context, the Boltzmann equation ([9]) is replaced by 

{v,Qx{fJ)) = (f I f{v)f{w)L{v')-ip{v)\dvdw). (19) 
\Jm.+ Jr+ l ^ I 

Here, we use the notation 

/oo 
h{s)d^[s). 

The unique possible collision invariants of the one-dimensional Boltzmann equation are 
obtained for '^{v) = 1 and ^{v) = v. 

The weak formulation of the Boltzmann equation can also be rewritten 

^{v)QxUJ)dv = \[ [ f{v)f{w){ip{v')+^{w')-^{v)-ip{w)) dvdw. (20) 

In (j20p the wealth variables u, w are nonnegative quantities, and the collision mechanism 
is given by (jl8p . A one-dimensional Boltzmann type equation of the form 

% = Qx{f,f) (21) 

based on the binary interaction given in (|18p has been considered in [18\ [23] and we 
refer to them for a deeper discussion of the model. Without loss of generality, we can 
fix the initial density /o(f) G 7^2(1^), with the normalization condition 



m{t) := / vf{v,t)dv = fh, (22) 
since by choosing <p{v) = v, (j20p shows that m{t) = m{0) for all t > 0. 
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As in section 12.21 one splits the collision mechanisms into a deterministic inelastic 
part and the random part : 

v' = V* + r]v; w' = w* + rfw 

where v* ,w* are deterministic wealth (corresponding to inelastic collision with constant 
restitution coefficient (1 — A)) 

y* = {I - X)y \W- W* =\V+{1-X)W. 

A formal Taylor expansion similar to ()2.2p . in the limit for A and r] small, leads 
to a drift term for the difference between {v,w) and {v*,w*) and a diffusion term 
proportional to the variance /?^. 

^{v') = ip{v*) + i]vdMv*) + ^V^v^d^ip"iv*) + ... 



Taking the average 



and on the other hand, the deterministic part gives 

ip{v*) = Lp{v) + \{w - v)dyip{v) + ^X^{v - wfdlip{v) + ... 

Inserting these expansions into the weak formulation of the Boltzmann equation p!9l 
and rescaling the time gives 



<{^,Qx{fJ))>= / f{v)f{w)\\{w-v)d,^{v) + \{\\v-wf+l3\^)dl^{v) 
12 L 2, 



dv dw. 



More precisely, the asymptotics of the one-dimensional Boltzmann equation for 
wealth distribution (|19p for A sufficiently small, and in the limit ^ — > 7, has been 
studied in [18] . In this so-called " continuous trading limit" , it is proved that the solution 
to the Boltzmann equation converges toward the solution to the Fokker-Planck equation 

which admits a unique stationary state of unit mass, given by the F-distribution 
where 

2 

, = l + ->l. 

This stationary distribution exhibits a Pareto power law tail for large velocities. We 
remark that in (j23|) the growth of the second moment follows the same law as the 
Boltzmann equation Q. 
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3 Quick overview of probability metrics 



In this section, we first briefly recall the main definitions and results about probability 
metrics and, more precisely, on Wasserstein {W2) and Fourier (dg) distances between 
two probability measures. 



3.1 Wasserstein distances 

Given two probability measures f,g ^ 'P(M^), the Euclidean Wasserstein Distance is 
defined as 

1/2 



W2{f,g) = inf \ [[ \v-x\^dUiv,x) 



(25) 



where H runs over the set of transference plans F, that is, the set of joint probability 
measures on x with marginals / and g G 'P(M^). From a probabilistic point of 
view, the Wasserstein distance can be alternatively defined as 

W2if,g)= mi _{E[\V-X\']f' (26) 
{v,x)er 

where F is the set of all possible couples of random variables {V, X) with / and g as 
respective laws. Let us remark that W2 is finite for any two probability measures with 
finite second moments f,g€z V2(^^)- 

The main properties of the Euclidean Wasserstein distance W2 are summarized in 
the following proposition. We refer to |13ll29l[3T] for the proofs and further information 
on the connections to optimal mass transport theory. 

Proposition 3.1 (VF2-properties) The space {V2i^^), W2) is a complete metric space. 
Moreover, the following properties of the distance W2 hold: 

i) Optimal transference plan: The infimum in the definition of the distance W2 
is achieved at a joint probability measure Hq called an optimal transference plan 
satisfying: 

Wi{f,g)= If \v-x\^dIio{v,x). 

ii) Convergence of measures: Given {fn}n>i md f in V2(^^), the following 
three assertions are equivalent: 

a) W2{fn,f) tends to as n goes to infinity. 

b) fn tends to f weakly-* as a measure and 

\v\'^ fn{v) dv ^ / \v\'^f{v)dv as n^+oo. 
Jm 

Hi) Convexity: Given f\, f2, gi and g2 in V2{^^) and a in [0, 1], 

W^iafi + (1 - a)f2,agi + (1 - 0)52) < aW^{fi,gi) + (1 - a)Wi{f2,g2). 
As a simple consequence, given f,g and h in V2(^'^), 

W2(,h*f,h*g)<W2{f,g) 
where * stands for the convolution in M^. 
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iv) Additivity with respect to convolution: Given fi, f2, gi and g2 in 'P-''™^^ 
with equal mean values, 

Wiih * f2,gi * 92) < <7i) + Wi{f2,g2). 

3.2 Fourier metrics 

Given / E P(M^), its Fourier transform or characteristic function is defined as 

/>)= / e-^--^df{v). 



Given any s > 0, tlie Fourier-based metric dg is defined as 

\f{k)-g{k)\ 

ds{f,g) = sup — (27) 

where = - {0}, for any pair of probabiHty measures f,g£ P(R^). This metric 
was introduced in [21] and further used in [16|. [T5l [28l I22j . Only recently, various 
applications to the large-time behavior of the dissipative Boltzmann equation [26^ [H [2] 
have revealed the importance of this distance. We refer to [T7] for a complete survey 
of this metric and the proofs of the statements below. 

The metric dg with s > is well-defined and finite for any two probability measures 
f,g Vs(^^) with equal moments up to [s] if s ^ N, or equal moments up to s — 1 
if s € N. The main properties of the ds metrics relevant to the ongoing discussion are 
summarized in the following result: 

Proposition 3.2 The distances dg with s > verify the following properties: 

i) Convexity: Given fi, f2, gi and g2 in VsO^^) with equal moments up to [s] if 
s ^ N, or equal moments up to s — 1 if s £ and a in [0, 1], 

ds{afi + (1 - a)f2,agi + (1 - 0)92) < a4(/i,ffi) + (1 - a)4(/2,52)- 

a) Superadditivity with respect to convolution: Given f\, f2, gi and g2 in 

Vs(^^) with equal moments up to [s] i/ s ^ N, or equal moments up to s — 1 if 

sen, 

dsifi * /2,9i *52) < ds{fi,gi) + ds{f2,g2)- 

4 Large time behavior for economy model 

4.1 Evolution of Wasserstein distance 

The Boltzmann equation (I2ip can be rewritten as 

df 

-Q^ = {fp+V * fq) - /> 
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where we use the shorthand fp{v) = {^/p)f{v/p) with p = X and q = 1 — X. Here, / 
is extended by to the whole of M in the convolution. The gain operator is defined as 
the measure given by 

{V>, QtU. /)) = (/ / •/'('^) /(^) (V^' ^{p+ri)v+qw) dv dw ) 

\Jr+ Jr+ I 

where is the Delta Dirac at the post-collisional velocity v' and (•,•) is the 

duality pair between continuous functions and probability measures. In probabilistic 
terms, the gain operator is defined as an expectation: 

Oa (/> /) =< /p+r/ * /g >= E \b{j,j^ri)V+qw\ 

where V and W are independent random variables with law / and independent with 
respect to the random variable ry. Here the expectation is taken with respect to all 
random variables. 

Let us take two independent pairs of random variables (y, X) and (VF, y) such that 
V and W have law /i and X and Y have law /2. From the convexity of VF| and the 
independence of the pairs, it follows that 

for any probability densities /i,/2 G 7^2 (IR)- Now, the last term is directly computed 
as the Euclidean distance between the two points {p + ri)V + qW and {p + ri)X + qY ^ 
and thus, 

^KOa (/i' h). Qt{f2, /2)) < E + 7?)(y - X) + q{w - y)|2] . 

Using independence of the pairs and taking the pairs to be optimal couples for the 
W2(/i,/2) in the probabilistic definition ([26l) . we deduce finally the property 

wHQtUi,fi),QtU2,f2))< [<{p + r]?>W] wKhJ2). 

Let us define, for s > 1 

&{s):={{p + r^r)+q' -I; (28) 

then 6(2) =< {p + rif > +q^ - I = 2A(A - 1) + /?^. It is not difficult to see that 
the convexity property of together with the Duhamel formula for (j2ip and the 
contractive estimate of the gain operator in W2 leads to the result: 

Theorem 4.1 Let fi{t) and /2(i) be two solutions of the one dimensional Boltzmann 
equation (j2T]) corresponding to initial values and in 7^2 (I^"^); satisfying conditions 
()22p . Then, for all times t >0, 

W2{fi{t)j2{t)) < exp{e(2)t}T^2(/?,/2°). (29) 

If (3'^ < 2A(1 — A), then S(2) < 0, and the Wasserstein metric decays exponentially to 
zero in time. 
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4.2 Evolution of Fourier metrics 



Analogous results for the evolution of the cis-metric (|27|) have been obtained recently 
in [23] by a suitable generalization of results in [25]. For the detailed computations we 
refer to [2^ . The study of the evolution of the metric (j27p , leading to the understanding 
of the large-time behavior of the solution to the kinetic equation (|2ip . requires a fine 
analysis of the quantity (j28p . As shown for the Wasserstein metric in the previous 
subsection, the sign of this quantity is in fact related to the contraction properties of 
the metric. Moreover, as has been noted in [24], the sign of (j28p is also related both 
to the number of moments of the solution which remain uniformly bounded in time, 
and to the possibility to conclude the existence and uniqueness of a steady state. The 
results in [2^ can be briefly summarized into the following 

Theorem 4.2 Take s > with &{s) < oo and let fi{t) and f2{t) he two solutions of 
the one dimensional Boltzmann equation ()2ip corresponding to initial values ff and 
in VrO^~^), satisfying conditions (j22p with r = max{l,s}. Then the following bound 
holds: 



where &{s) is given by (f28]) . 

Also, the temporal behavior of the moments is almost completely determined by 
the function &{s). 

Theorem 4.3 Let s > 1 and /o G VsO^~^) with < ©(s) < oo and let us denote 



Then, for the weak solution to the Boltzmann equation, the following estimates hold: 
1. If ©(s) > 0, then, as t ^ oo, 



2. If &{s) < 0, then the sth moment is bounded for all times. Moreover, as t ^ oo. 



Here, the remainder terms o(l) converge to zero exponentially fast. 

Another important conclusion of the analysis of [21] is that the essential function 
©(s) does not only decide whether or not the steady state foe develops a Pareto tail. 
In fact, the positive zero of ©(s) actually determines the value of the Pareto index. 

A comparison of the contraction results for the Boltzmann equation (I2ip shows 
that the contraction properties are heavily linked, through the key function (|28p . to 
the (eventual) formation of tails. While the situation for equation (|2ip is reasonably 
well understood, the corresponding analysis for the Boltzmann equation ([9]) deserves 
further investigation. We will discuss equation Q in detail in the following section. 



ds{fi{t),f2{t)) < exp{e{s)t} d,{flf^), 



(30) 
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5 Large time behavior for stochastic granular media 

Let us consider here the modification of the Inelastic Maxwell Model introduced in [6] 

% = QeifJ), (31) 

where the collision operator is defined weakly as 

i^,Qe{f,f)) = ^( [ I I fiv)fiw)Liv')-^iv)\dadvdw). (32) 
4vr \Jk3 Jj^3 Js2 I J / 

As discussed in the introduction, the collision mechanism relies on a random coefficient 
of restitution, 

/ 1, , 1-e 1 + e,, 
^ = 9(^ + '"^)"l 

(33) 

' 1, , 1-e l + e,, 
w = -{v + wj : — u : — \u\a. 

As before, we write u = v — w, e = e + ij and r/ is a real-valued random variable, with 
zero mean and variance /3^, given by a measure /i(s) with support on [— e,oo). Here, 
< • > means the expectation with respect to i.e., the integral over M with respect to 
fi. 

It is quite straightforward to check that conservation of mass and momentum re- 
mains and that 

/I '|2 I I /|2 I |2 I |2\ n 

[\v \ + \w \ — \v\ — Iw] ) = 

for the model that is conservative in the mean in which = 1 — e^. From (|32p . we 
deduce that the temperature evolution is 



^ / \v\^f{t,v)dv = 0, 
at J^a 



and thus we deduce that 6{t) = 9{0) for all times t > and we will fix it to one for 
convenience. 

5.1 Evolution of Wasserstein distance 

Given a probability measure / on M'^, the gain operator is in fact a probability measure 
QtifJ) defined by 

{'P,Qt{fJ)) = ( [ I f{v)f{w){ipMv,u,,r,)dvdw 

where Uv^w,ri is the uniform probability distribution on the sphere S^^w with center 
Cv,w = ^{v + w) + [^-^ + r]j [v — w) and radius r^^^ = ^-^\v — w\ as in [H]. In 
probabilistic terms, the gain operator is defined as an expectation: 

where V and W are independent random variables with law / and independent of the 
law of rj. As in [H], we get the following result: 
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Theorem 5.1 Given f and g in 7^2 (IR^) with equal mean velocity, then 

W2{Qt{fJ),Qti9,9))<W2{f,g). 

Proof.- Let us take two independent pairs of random variables {V,X) and (M^, y) 
such that V and W have law / and X and Y have law g. Also, let us take two 
independent random variables r] and fj with law fi. Convexity of W2 implies 

WiiQtif, f),Qti9, 9)) = WiiE [Uv,wA , E [Uxxd) < E [Wi{Uv,w,v^Ux,Y,v)] (34) 

where the expectation is taken with respect to the joint probability density in M^^ of 
the six random variables. Here, the independence of the pairs of random variables has 
been used. 

As proved in [13], the VF| distance between the uniform distributions on the sphere 
with center O and radius r, Uo,r, and on the sphere with center O' and radius r', Uo',r', 
in is bounded by \0' - Op + (r' - rf. 

We now estimate the right-hand side of (j34p by using the formulas for the center 
and radii of the spheres given in (j33p to deduce 

W^iQtU, /), Qt{9, 9)) < ( ^"'g^'' ) E[\V-X\^]+ (^^) K[\W-Y\'] 

+ (L^^E[iV-X).{W-Y)] 

where the Cauchy-Schwartz inequality has been used. 

Finally, we take both pairs {V, X) and {W, Y) as independent pairs of variables with 
each of them being an optimal couple for VF2 (/)<?) in the probabilistic definition (j26p 
to obtain 

WHQtif, /), Qt{9, 9)) < ^±^±J^ WIU, g) + E [(y - X) ■ {W - Y)] , 

where the last term is zero because the random variables are independent and have 
equal means. Since = 1 — in the conservative case, the result is proved. □ 

As a consequence of the previous property of the gain operator, we draw the fol- 
lowing conclusion about controlling the distance between any two solutions of ([3T]) in 
the conservative case. 

Theorem 5.2 If fi and f2 are two solutions to (|3ip with respective initial data /{* and 
/2° in 7^2 (IR^) with zero mean velocity, then, for all t > 0, 

wKh{t)J2{t))<wl{flf^2). 

Proof.- Duhamel's formula for (j3ip reads as 

/^(^) = e-*/°+ fe~^'~''>Qt{fi{s),f^{s))ds, i = l,2. 
Jo 
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As before, the convexity of the squared Wasserstein distance in Proposition 13.11 and the 
contraction of the gain operator in Theorem 15.11 imply 



Wi{h{t),f2{t))<e-'wi{f'^J^)+ / e-^'-'^Wi{Qt{fi{s),h{s)),Qt{f2{s),f2{s)))ds 

Jo 

< e-* W^K/O, /O) + / e-(*-^) Wiihis), Ms)) ds. 
Jo 

Therefore, the function y{t) = e* VFf /2(i)) satisfies the inequahty 



y{t)<y{0)+ / y{s)ds 
Jo 

and thus y{t) < y(0) e* by Gronwall's lemma, concluding the argument. □ 
5.2 Evolution of Fourier metrics 

We start by writing a closed form of the Boltzmann equation in Fourier variables. In 
fact, it is not difficult using Bobylev's identity in O [U O [6] to get 

1 



where 



Qtif, f) = ^{ fit, k.)fit, K) da 



1 + e l + e,,, , , 3-e l + e 

-k ■ — \k\a and kj^ = — — k -\ ; — \K\(J . 



4 4 ' ' '4 4 

Let us start by analyzing the evolution of the distance ^2 that in view of the properties in 
Propositions 13. l l and 13.21 should verify the same non-strict contraction as the transport 
distance W2. 

Theorem 5.3 Given f and g in 7^2 (1^^) with equal mean velocity, 
d2iQtifJ),Qti9,9)) < ^-^^-^d2{f,g). 



Proof. - Using the Fourier representation formula above, we deduce 

' f{k^)f{k^)-g(k^)g{k+) 



Qt{f,f){k)-Qt{g,9){k) _ 1 

" Air \ .Ig2 



for all k G Mo- We now estimate the integrand as 

f{k.)f{k+)-g{k.)g{k+) 



da 



<SUD|M^lM\rM±^ 



d2if,g) 



and thus 



d2{QtifJ),Qti9,9))< 



52 



dcr ) d2if,g). 
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We observe that 



is a function of the angle between the unit vectors k/\k\ and a and the random variable 
r/, and that 



^ ' 47r 



■ da 



52 



3 + + 



In fact, we can compute 



1 + e 



2 1 - cos 1? 



3 — e 



+ 



1 I 2 
1 + 



+ 2 



3-e\ /1 + e 



cos t9 



(35) 



where ■!? is the angle between the unit vectors k/\k\ and a from which the value of / 
is obtained. Putting together previous estimates we get the contraction in d2 with the 
same constant as VF| as desired. □ 



Now, let us see that we can also control Fourier-based distances with exponent 2+ a, 
a G [0, cxd). Let us set 



2t(a,e,r?) := 



1 



+ 



1 I ~\ 2 
1 + 



2(1 - COST?) 



2 + a 
2 



3-e 



+ 



1 + e 



+ 2 



3-e\ /1 + e 



cos'd 



2+a 
2 



sin •§ d-d 



4 + a 



l + e\2+"^l_ 



l_e|4+a' 
2 I 



1 



2 I 



(36) 



Whenever there is no confusion, i.e. for e and rj fixed, we will denote just by 2t(a) the 
above constant. 

Theorem 5.4 Given f,g£ 'P2+a(M^) with equal moments up to order 2 + [a], there 
exists an explicit constant 2l(a, e, 77) > given by ([36]) such that 

d2+aiQt{f, f),Qt{9, g)) < 2i(«, e, r,) d^+M, g)- 
Proof.- As in the proof of the previous theorem, we compute 



Qt{f,f){k)-Qt{g,g){k) 



\k\ 



2+a 



1 



f{k+)f{k-)-g{k+)g{k-) 



< A sup 

fcGK§ 



52 \k\ 

\fik)-g{k)\ 



2+a 



da 



\k\ 



2+a 
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where A is given by 



■^^=sU. — w^' — "')■ 

By inserting the expressions of k_ and into (j37p and computing the integral we 
conclude ^ = 21(q, e, ry) and the proof follows. □ 



As a consequence, we obtain an estimate on contraction/expansion of the Fourier 
distances d2+a between solutions. 

Theorem 5.5 Let a > be such that 2l(a,e,r/) < oo. Let fi and /2 be two solutions 
to (|3ip corresponding to initial values fi, f2 with equal moments up to 2 + [a]. Then, 
for all t > 0, 

d2+Ml{t)j2{t)) < fi2+a(/{',/2°)e-^("'^''')*, (38) 

with C{a, e,r]) = 1 — 2l(a, e,r]). 

Proof.- The Fourier expression of equation ()3ip is given by 

% = ^ l\ fik+)hk.)da -f = QtiLf) - f, 

whose solution satisfies 

/(t, k) = e-*/(0, k) + /* e-(*-^)Q+(/r/)(s, k) ds . (39) 
Jo 

Taking the expressions of the two solutions fi{t) and /2(i) in (I39p . subtracting them 
and dividing by jfcp"'"" with /c G Rq, we get 

t Ul-h)it,k) _ fl{0,k)-h{0,k) f ^ (Q^(/l,/l)-Qe+(/2,/2))(^,fc) 

Using Theorem 15.41 and taking the supremum in k £ Mq, we obtain 

e'd2+a{Lhm < d2+a(/i(0),/2(0)) +2t(a,e,r/) /" e^d2+a(/i, /2)(s)ds. 

JO 

Let us set w{t) = e*(i2+a(/i, /2)(i)- Then 

w{t) < w{0) +^{a,e,r]) [ w{s)ds, 

Jo 

which by Gronwall's inequality implies w{t) < ii;(0) e^*^"'*^'^)*, concluding the proof. □ 

The function 2t(a) : [0, oo) — > is convex by direct inspection. Taking into ac- 
count that 21(0) = 1, there are only three possible scenarios for the qualitative behavior 
of 21. These are characterized by the sign of 21' (0). In case 21' (0) > 0, the function 21(a) 
has a minimum at a = due to convexity, and thus 2t(a) > 1 for all a > 0. In this 
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case, there does not exist any a G M+ such that 2t(a) < 1 and there are no contraction, 
only expansion, estimates of dg for s > 2. 

Suppose that 2l'(0) < 0. In this case, the contraction properties of ds depend on 
whether 

hm 21(q) < 1 

or 

hm 21(a) > 1. 

In the former case, 2t(a) < 1 for q > 0. Theorem 15.41 then imphes that the d^-metric 
is contractive for all values of the parameter s > 2. In the latter, since 21(0) = 0, the 
convex function 21(a) has a minimum attained at some point a > 0, and at the same 
time there exists a > a for which 21(a) = 1. Thus, 21(a) < 1 in the interval < a < a, 
and at the same time 21(a) > 1 for a > a. In this case Theorem 15.41 implies that the 
Boltzmann equation is contractive up to but not including order a. 

Remark 5.6 In order to clarify the behavior of^{a,e,rj), we can fix the random vari- 
able rj to assume only two values, while respecting conditions ([5]). This can be done 
by assuming that rj only takes the value \/l — e'^/g with probability g'^/{l + g"^) and the 
value Vl — e'^g with probability 1/(1 + g"^). By varying the parameters g and e one 
encounters the whole variety of possible behaviors of the function 2l(a,e,r/). Since 



2t(a,e,7?) = ^— ^ + 



2 



4 + a\ V2/ i_ \l=e\^ / 

2t(a,e,r/) results in the sum of four contributions, one of which is 

^ ) ■ 

For any fixed values of a > and e, since the numerator grows like g'^~^°', we can choose 
g » 1 in such a way that C{a,e,ri) > 1, and Theorem 15.41 implies that the Boltzmann 
equation is contractive up to but not including order a. 

On the other hand, choosing for example a = 2 to simplify computations, one 
obtains easily 

. + fl-eV fl-eY\ 23 - e + (e^) 



Choosing now 1 — e << 1, and g = \/l — e^/e, one obtains that e assumes the value 
with probability 1 — and the value 1/e with probability e^. Therefore (e^) = 1/e^, 
which implies 2t(2,e,r/) < 1 as long as 1/e^ — e < 2. In this second case Theorem 15.41 
implies that the Boltzmann equation is contractive at least up to order 4. 
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5.3 Existence and uniqueness of regular isotropic steady states 

Existence and uniqueness of steady states, as well as the size of their overpopulated tails, 
can be derived in full generality (that is, without imposing restrictive conditions on the 
random coefficient of restitution) by adapting to the present situation the methodology 
of [2], which refers to the inelastic Boltzmann equation for Maxwell molecules. This 
methodology, in fact, is based only on the contractivity properties of the d^-metric, 
which are analogous to Theorems 15.41 and 15. 5[ 

It has to be remarked that the approach in [2] is not suitable to recover the (eventual) 
regularity of the steady profile. A regularity result for the steady state of the inelastic 
Boltzmann equation for Maxwell molecules has been obtained in a recent paper by 
Bobylev and Cercignani In this paper they were concerned with properties of the 
self-similar profiles of the Boltzmann equation for both elastic and inelastic collisions, 
and, in addition to the existence, they obtained results on the regularity of the steady 
profiles by showing that the Fourier transform of the steady profile satisfies a suitable 
upper bound. Their method takes advantage of the existence of a super-solution to 
the rescaled equation in Fourier variables (BKW-mode). In our collisional setting, the 
situation is more involved, and it requires a precise analysis. 

In Fourier variables, the steady state of ([9]) is a solution of the integral equation 

1 



f{k_)f{k+)da)=f{k), (41) 

52 



where fc+ and k- are given by the relations 

k- = — ^^/c ^^|A;|(T and k^. = - — - k H \k\a . 

4 4 4 4 

Since isotropy is not destroyed by the collision operator, by choosing isotropic initial 
values, one concludes with the isotropy of the (eventual) steady state. Taking this 
property into account, the following result can be obtained as a consequence of Theorem 
15.51 (see j2] for details) . 



Corollary 5.7 Equation (j3ip has a unique isotropic steady state foo in the set of 
isotropic probability measures with unit mass, zero mean velocity and unit temperature. 
Moreover, given any solution f to ([3T|) for the initial data fo G V20^^) with zero mean 
velocity and unit pressure tensor, 



d2+a{f{t), foo) < d2+aifo, foc)e 



-C{a,e,T])t 



for all t > 0, < a < 1. Thus, if^{a,e,r]) < 1, f{t) converges to the stationary state 
as t ^ 00 in the d2+a sense. 

Remark 5.8 The previous result shows that the stationary states attract all solutions 
with initial data having zero mean velocity and unit pressure tensor. The assumption 
of having unit pressure tensor can be weakened to having initial unit temperature by 
proceeding similarly to the homogeneous cooling state analysis in [HI [2] . 
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Let us define 



b\e,r^,e) 



(42) 



Recalling the definition of A;_|_ and A;_ given in (j35p . it is immediate to show that 

a + b>l; ]^l^j^{a^ + h^) sinOdd^ = 1 (43) 

The first property in ()43p is a direct consequence of the equality kj^ + = k, while 
the second is the equality 21(0) = 1 in (j37p . Let us set x = \k\. Then, for any function 
il^{x), the left-hand side of (14ip can be rewritten in the form 

R[il;{x)] = \(^j^ ij{ax)ip{bx) smOde"^ . (44) 

Under the conditions of Corollary 15.71 the Boltzmann equation has a unique steady 
state foo{x), of unit mass, zero mean velocity and unit second moment. 

Let us remark that < R[tp] < 1 if <'(/'< 1, and Rlip] < R[(j)\ if < V < 0- 
Hence the iteration is monotone increasing and converges point-wise if we choose the 
initial approximation < (/^o < 1 in such a way that <j)Q < R[(j)o]- As observed in [8], 
4>o{x) = exp{— 2:^/2} allows us to obtain a monotone increasing sequence. In fact, since 
the function e~'^ , r > is convex, by Jensen's inequality we obtain 

This implies that the limit foo{x) > 0. The trivial limit /oo(a;) = 1 will be excluded if 
there exists a non-zero function (j)o{x) such that 

M^) < V'o(^), (46) 

and at the same time ipo{x) generates a monotone decreasing sequence. 

Inspired by the ideas of Desvillettes et al in [19], given a fixed positive constant p, 
we introduce the fixed point operator 



foo{x) llx < p 

R[4j{x)] ifx>/9 



on bounded complex functions : M — > C. Notice that R is closely related to the 
Fourier transform of the collision kernel. 

Lemma 5.9 Let < foo{x) <1 he the steady state of the Boltzmann equation, and let 
us define 



foo{x) if x < p 
exp(— ^x) if X > p. 
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Then, if the random variables a{e, r], 9) and b{e, rj, 6) are such that 

P{a < 6) + P{b < 6) ^ as (5 ^ 0, (47) 
there exist positive constants p and fi such that 

Proof.- Since the steady state is an isotropic probability density function of unit 
mass, zero mean velocity and unit second moment, there exist positive constants M < 
1/2 and p such that (cfr. [19]) 

< foc{x) < e"*^^' if x<p. (48) 

Hence, we can fix p and M to obtain 

Mx) < e-^""^ if x<p. (49) 

Clearly, thanks to the definition of if ^ there is nothing to prove. Therefore, 
let us consider the possible cases corresponding to x > p. Since a + 6 > 1, if both 

ax > p, bx > p, 

{i,o{ax)Mbx)en < 1- 
If now both ax < p and bx < p, using bound (jl9]) . we obtain 

{M<^x)^o{bx)en < (es(^)), 

where 

g{x) = px- M{a^ + 6^)x^. 
Since a + 6 > 1, it follows that a^ + b"^ > 1/2. Thus 

g{x) < px - ^{a'^ + b'^)x'^ < if /i < (50) 

Consider now the case in which ax < p, while bx > p. In this case 

where 

h{x) = p{l - a)x - Mb^x"^. 
Since a + 6 > 1, it follows that b> 1 — a, and 

u? 

h{x) < z{bx) = pbx - M{bxf < (51) 

In fact, the function z{r) has a maximum at f = p/[2M). Moreover, since z{r) decreases 
for r > f, if r > 3f, 

z(r)<z(3r-) = -3^. (52) 
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Let us split the calculation of the mean value into the sets A = {bx > 3r} and A'^ = 
{bx < 3f}. Thanks to conditions ([5T|) and ([52]) one obtains 

(e'^(^)) < P{A) exp {-3^} + P{A^) exp {^} • (53) 

Let us set S = 3f = 3/i/(2M). By hypothesis, since x > p, 

P(A^) = P{bx <6)< P{bp < 5) ^ if (5 ^ 0. 
Consider that we can rewrite ([53|) as 

(gM^)) < (l-P(^^)exp|-^M,52| +P(^^)exp|^M,52 

1 - -(1 - 4P(A'=))M(52 + o(j2) < 1 
3 

if 6 is sufficiently small. Now, this condition on 6 can be satisfied by choosing p 
sufficiently small. This is not in contrast with condition ()50p . Since the case in which 
ax > p while bx < p can be treated likewise, the lemma is proven. □ 

Remark 5.10 Condition ()47p excludes some pathological situations related to the def- 
inition of the random variable ij that describes the randomness of the coefficient of 
restitution e. For example, condition (|47p is violated if i] is concentrated on some 
particular point, 

P(n = 1 - e) = p > 0. 
In this case, in fact, P{b{e, rj, 0) = 0) = p, and condition (|47p is false. 

Lemma [5.91 implies that, starting from ipQ, the iteration process leads to a monotone 
decreasing sequence. On the other hand, it is clear that, for p sufficiently small, 

< M^) < M^) < 1- 

Given ^ > 0, define as the set of functions ip with -0(0) = 1, ^'(0) = /^dCO), and 
satisfying the estimates 

< exp(— Kx^) for X < p, \'^{^)\ < exp(— ^x) for x > p. (54) 

The previous inequalities prove the following 

Theorem 5.11 For any pair of functions a and b satisfying conditions (j43p anti (j47p . 
the integral equation (|44p has a nontrivial solution fooi^) such that foo{^) belongs to 
the Gevrey class K^^ defined by 

Remark 5.12 An analogous regularity result can be proven for the steady state to the 
one- dimensional kinetic model (|2ip 124^ . In this case, it is important to know that the 
mean wealth of the stationary state is equal to one. 
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5.4 Fat tails of stationary states 

In this work, we will only examine the case of the fourth moment, postponing the 
complete analysis of moment evolution to future research. Here, we will show that 
under certain conditions on the random variable, the fourth moment diverges or is 
controlled uniformly. 

Lemma 5.13 Let the restitution coefficient e and the random variable i] be chosen so 
that 21(2, e,7y) < 1. If f^ is a Borel probability measure on M'^ such that 

I \v\'^ f^{v) dv < oo, 

then the solution f to ([3T]) with initial datum f^ satisfies 

sup / {vl"^ f{t,v) dv < oo. 
r>0 Jr3 

Proof.- Without loss of generality we can assume that f'^, and hence f{t) for all 
r > 0, has zero mean velocity and unit temperature. We let 

m4{t) = / Ivl"^ f{t,v) dv 

denote the fourth order moment of f{t). Then, using the weak formulation of the 
inelastic Boltzmann equation, we have: 



dmi{t) 



v\^Qe{f{t)J{t)){v)dv (55) 



dt 

that can be computed as in [T3| by 

\v\^Qe{f,f)iv)dv=-<C> [ \v\^fiv)dv+<f,,> ( [ \v\^f{v)dv) 

Jr3 ^ Jr» ' 

+ < /X2 > // {v ■ wf' f {v) f (w) dv dw 



2 



where 

1 1 

/"I = ^('^i + - J^s) and [i2 = -{i^i - V2) 



with 



Z., = (e2 + e'2)2_l + leV2, z.2 = 2[62 + 6'2-l + V], ^3=4(6^-1), 



and 



= 1(1 + 4e - Te^ + 46^ - 2e'^) with e = - — - and e' = 1 - e. 
3 2 



Now, ([55]) reads 

^^^^ = - < C > m^it) + m{t) (56) 
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where m{t) is a combination of second order moments, which are bounded in time 
since the kinetic energy is preserved by equation (j3ip . Moreover one can check from 
the expression of C in terms of e that < ( >= 1 — 21(2, e,ry) > 0. This ensures that 
m4(t) is bounded uniformly in time if initially finite, and concludes the argument. □ 

The preceding result also shows the divergence of the fourth moment in case the 
random variable r] and the restitution coefficient e are chosen to satisfy 21(2, e,rj) > 1 
but 2t(Q, e,ri) < 1 for some < a < 2. 

Corollary 5.14 Let the restitution coefficient e and the random variable r] be chosen 
so that 21(2, e, r/) > 1 but 2l(a, e,rj) < 1 for some < a < 2. Then, the unique isotropic 
steady state f^o in 7^2 (IR'^) of equation (|3ip with zero mean velocity and unit pressure 
tensor has unbounded fourth moment. 

Proof.- With the notation of the previous subsection, the evolution of the fourth 
moment for isotropic densities given in Lemma l5 . 1 3 1 ensures that 

= - < C > m4(0 + m(t), 

where m[t), which is a combination of second order moments, is bounded from below. 
Recall that < C >= 1 - 2t(2, e, ??) < to conclude. □ 
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